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ABSTRACT 

We analyze the star formation history (SFH) of the Galactic disk by using an in- 
fall model. Based on the observed SFH of the Galactic disk, we first determine the 
timescales of the gas infall into the Galactic disk (tin) and that of the gas consumption 
to form stars (tsf). Since each of the two timescales does not prove to be determined in- 
dependently from the SFH, we first fix tgf . Then, tin is determined so that we minimize 
X^- Consequently, we choose three parameter sets: (tgf [Gyr], tin [Gyr]) = (6.0, 23), (11, 
12), and (15, 9.0), where we set the Galactic age as 15 Gyr. All of the three cases predict 
almost identical star formation history. Next, we test the intermittency (or variabil- 
ity) of the star formation rate (SFR) along with the smooth SFH suggested from the 
infall model. The large value of the statistic supports the violent time-variation of 
the SFH. If we interpret the observed SFH with smooth and variable components, the 
amplitude of the variable component is comparable to the smooth component. Thus, 
intermittent SFH of the Galactic disk is strongly suggested. We also examined the 
metallicity distribution of G-dwarfs. We found that the true parameter set lies between 
(tgf [Gyr], tin [Gyr]) = (6,23) and (11, 12), though we should need a more sophisticated 
model including the process of metal enrichment within the Galactic halo. 

Subject headings: Galaxy: evolution — Galaxy: stellar content — methods: statistical 
— stars: abundances — stars: formation 



1. INTRODUCTION 

In the field of the Galactic chemical evolution, Eggen, Lynden-Bell, & Sandage (1962) have 
inspired the modeling. From the correlation between the ultraviolet excess and orbital eccentricity 
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of stars, they concluded that the Galaxy formed by collapse on a free-fall timescale from a single 
protogalactic cloud. An alternative picture of halo formation is proposed by Scarlc & Zinn (1978). 
They argued that Galactic system is formed from the capture of fragments such as dwarf galaxies 
over a longer timescale than that proposed by Eggen et al. In any case, determining the timescale of 
the infall of matter and the chemical enrichment is an important problem to resolve the formation 
mechanism of the Galactic system. Indeed, there have been a number of papers that investigated 
the formation and chemical evolution of the Galaxy (e.g., Matteucci & Prangois 1989) and spiral 
galaxies (e.g., Lynden-Bell 1975; Sommer-Larsen 1996). 

Many models of the star formation history (SFH) of the Galaxy includes the formation of the 
Galactic disk by gas infall. This scenario (so-called the infall model) is consistent with the age- 
metallicity relation of the disk stars (e.g., Twarog 1980), if a reasonable SFH is applied. Moreover, 
the infall model provides a physically reasonable way of solving the G-dwarf problem (e.g., Pagel 
1997, p. 236), contrary to the closed-box model which tends to overpredict the number of the low- 
metallicity stars. 

One of the factors that determine the star formation rate (SFR) is the gas content of galaxies. 
Indeed, the SFR and the gas density is closely related to each other (Kennicutt 1998), and this 
relation is generally called the Schmidt law (Schmidt 1959). The SFR and the gas density, p, are 
related as SFR cx /o", where n = 1 — 2. In the one-zone model, the relation is simply assumed to 
be described with SFR oc Mg , where Mg is the gas mass in the galaxy. 

Though the "classical" infall model is widely accepted, there are observational data that suggest 
the intermittent star formation activities in spiral galaxies. Kennicutt, Tamblyn, & Congdon (1994) 
showed the ratio of present-to-past SFR in spiral sample has a significant scatter. More recently, 
Tomita, Tomita & Saito (1996) analyzed the far-infrared to B-band flux ratio /fir//b of 1681 
spiral galaxies (see also Devereux & Hameed 1997). The indicator /fir//b also represents the 
ratio between the present SFR and the averaged SFR during the recent Gyr. They showed the 
order-of- magnitude spread of /fir//b and suggested the violent temporal variation of the SFR. 

The intermittency of the SFH in the Galactic disk is recently suggested by Rocha-Pinto et al. 
(2000; hereafter ROO). They derived the SFH of the Galactic disk from observed age distribution 
of the late-type stars and suggested that the star formation activity of the disk is intermittent or 
violently variable. They tested to show that the variation is of statistical significance. Their test 
are based on the null hypothesis that the SFR of the Galaxy is temporally "constant", but their 
derived SFH indicates that the SFR around 10-15 Gyr ago is significantly smaller than that around 
0-10 Gyr ago. This trend of the Galactic SFR can be interpreted along with the infall scenario of 
the Galactic disk formation: Gas-infall to the Galactic disk occurred in the first ~ 5 Gyr. Thus it 
is an important work to interpret their data through the infall model and to examine whether the 
variation of SFR is still significant even in the infall scenario. 

In this paper, we re-analyze the data by ROO in the context of the infall model. In the next 
section, we briefiy review the derivation of the SFH in the Galactic disk by ROO. Next, in §3, we 
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present the formulation of the infah model, based on which we test the data by ROO statistically. 

In §4, we interpret the observed SFH by two components; underlying infall-model component and 
violently variable component. In the same section, possible physical mechanisms for the variation 
of star formation are also presented. Finally, we summarize the content in §5. 

2. THE DATA 

ROO provided the SFH of the Galaxy inferred from the stellar age of the solar neighborhood. 
They used 552 late-type stars. The age of each star is estimated from the chromospheric emission 
in the Call H and K lines (Soderblom, Duncan, &; Johnson 1991). After metallicity-dependent age 
correction, completeness correction, and scale-height correction^, they derived the age distribution 
of the stars. Then, after correcting the stars that is not alive at present, they derived the SFH of 
the Galactic disk in their Figure 2. The discussion in this paper is based on Figure 2 of ROO. 

ROO concluded that the inferred SFH is representative of the SFH of the whole disk, since the 
timescale of the diffusion of stars to the kpc scale is ~ 0.2 Gyr, which is much shorter than the 
Galactic age. Thus, in the next section, we model the Galactic SFH by a one-zone model to extract 
the global property of the SFH. The difference between the one-zone and the multi-zone treatment 
is seen in Figure 6 of Sommer-Larsen (1996), where we see that the one-zone treatment is a good 
approximation. 

3. SFH FROM INFALL MODEL 

In this section, we give a physical interpretation of the SFH by ROO. First, we adopt the infall 

model for the interpretation, since the model has been successful in reproducing the age-metallicity 
relation of stars in the solar neighborhood (e.g., Pagel 1997). Next, we infer the timescales of 
the following two processes; the infall of gas into the Galactic disk and the gas consumption to 
form stars. Finally, we examine whether the intermittent SFH proposed in ROO is of statistical 
significance. 

3.1. Model Description 

We assume the one-zone model: The gas mass in the Galaxy is treated as a function of time, 
Mg(t). The time evolution of the gas mass is determined by the infall from the halo, whose rate 
is described by F{t), the consumption by star formation, ip{t), and the recycling from stellar mass 
loss. If the instantaneous recycling approximation is adopted (Tinsley 1980), the time evolution of 
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the gas mass in the Galaxy is described by 



'^^^'--{1-R)^ + F, (1) 



dt 

where R is the returned fraction from stellar mass loss, described by 



R = {m — Wm)<p{'>Ti) dm . (2) 



mt 

Here, ^(m) is the initial mass function (IMF) of stars, mt is the present turnoff mass (IMq), 
m„ is the upper cutoff of the stellar mass (lOOM©), and Wm is the remnant mass (we assume 
Wm = O.7M0 for m < 4Mq and Wm = IAMq for m > 4:Mq). The IMF is normalized so that the 
integral of m(j){m) in the full range of the stellar mass becomes unity. When the Salpeter IMF, 
0(m) oc m~^'^^ (Salpeter 1955) with the lower cutoff of 0.1 and the upper cutoff of 100 Mq is 
assumed, R = 0.32. Pagel (1997) derived a similar value for R by using a different form of IMF 
and remnant mass {R = 0.2 - 0.3), and commented that the uncertainty in i? is ~ 0.1. 

Since the normalization of the SFR and statistical test become much more complicated unless 
the SFR is presented by analytic function, we adopt the instantaneous recycling approximation for 
clarity in this section. The analytic form makes it significantly easy to perform the statistical test 
in §3.3. We examine the propriety of adopting the instantaneous recycling approximation to this 
problem in §4.2. 

In this paper, the timescale of the gas infall onto the Galactic disk, tin, will be determined by 
a fitting to the observational data in §3.3. For the convenience of the fitting, we assume the infall 
rate is expressed by an analytic function. A natural form is the following exponential function 
(Pagel 1997, p.242): 

F{t) = ^eM-t/tin), (3) 
where Mq indicates the total mass that can fall into disk; in other words, 

I 

F{t)dt = Mo. (4) 



/ 
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Normalizing the equation (1) by Mq leads 



^ = -{l-R)k + l. exp (-iAin) , (5) 

a t tgi tin 



where 



The Schmidt law of n = 1 (§1) is assumed as follows: 

^ = ^/Mo = /g/trf , (7) 
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where tgf indicates the timescale for the gas to be converted to stars (i.e., the gas consumption 
timescale). Then, we can express the solution of equation (5) analytically as 

= j^{e-"«--e-<'-'»"'-) , (8) 

One of the goals in this paper is to infer the two timescales, iin and igf from the SFH in ROO. Thus, 
we define the SFR at the Galactic age of t, ip, as 

lij = i^(t;tin,tsi) = fg/tsi- (10) 

For the following discussions it is convenient to define the averaged SFR for the normalization 
of i/j. The averaged SFR, (t/;), is defined by 

(^) = ^_^^°Vidt. (11) 
Equations (8) and (10) lead the averaged SFR as 

(^) = ± r° ^ L-tlh^ _ e-(i-«)*/tsf ) di, (12) 



where Tq is the age of the Galactic disk and we take Tg = 15 Gyr following ROO. Then we obtain 

^ g-T/Tin _ g-(l-R)T/rsf 



(V') Tin (1 - e-Vrin) _ (1 _ e-(l-«)/r.f) 

where r = t/To, Tin = tin/Tb, and Tgf = tgf/^G- Here we note that 



= *(r; r,„,rsf) , (13) 



r*(r) 

JO 



dr = 1 . (14) 



3.2. Trend Estimation 

We extract the overall "trend" from the data of ROO. We used the smoothing method developed 
in the field of the exploratory data analysis (EDA), as well as the ordinary moving average. The 
EDA smoothing is based on the moving median, which is is known to be quite robust against outliers 
in the datasets compared with the moving average (e.g., Hoaglin, Hosteller, &: Tukey 1983). Since 
the SFH in ROO violently varies with time, the EDA procedure is expected to be suitable for 
the problem. Detailed procedures are summarized in the appendix. The smoothed results with 
the above two procedures are shown in Figure 1. The original data are depicted by the dotted 
histogram. In the upper panel, we show the smoothing results by moving average. Details of the 
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moving average method is extensively discussed in e.g., Hart (1997). The dot-dot-dot-dashed hne 
is the smoothed SFR with smoothing kernel width of 1.2 Gyr, strong dashed line is the smoothed 
SFR with 2.0 Gyr kernel, and dot-dashed line is the smoothed SFR with 2.8 Gyr kernel. In the 
lower panel, the dashed line represents the moving average of 2.0 Gyr kernel, and the solid line is 
the EDA smoothed SFH. We observe that the 1.2- Gyr kernel is not sufficient to smooth the varying 
SFH. This means that the timescale of the variation is less than ~ 2 Gyr. Possible mechanisms of 
the variation that works in less than 2 Gyr are listed in §4.4. 

We see that the two results yield consistent trends, but the EDA result is not affected by the 
sudden jump of the values compared with the moving average as clearly seen at Tq — t = 10 - 
15 Gyr. We find the rise of the SFR in the early epoch. This is interpreted as being the gradual 
increase of gas by the infall. Thus we suggest that the overall SFH of the Galaxy is well described 
by the infall scenario. We statistically infer the infall timescale and the gas consumption timescale 
in the next section. We will also examine whether the 'residual' deviated from the trend is of 
statistical significance there. 



3.3. Parameter Estimation and Statistical Test 

In this subsection we perform a statistical test for the observed Galactic (normalized) SFR, 
^'obs (shown as SFR/(SFR) in ROO) with respect to that inferred from the infall model. The 
observed data are binned as ^'obsj) i = 1,2, . . . ,k. Once the parameter set (rjn, Tgf) is fixed, we can 
estimate these parameters from the modified Pearson's chi-square statistic, xo^ (see e.g., Rao 1973, 
p.352): 



where is the raw number of stars in the z-th bin, N is the sample size (in this case N = 552), ^'Jj^^j 
is the raw normalized SFR without completeness corrections, and is the binned theoretical 
SFR which we would have observed under the same condition as ^'^fj^j- However, since the SFR in 
Figure 2 of Rocha-Pinto et al. (2000a) is corrected for the selection biases when they have derived 
the value, we must include the effect of the incompleteness correction in the statistical analysis. 
We set ncor,i = CiUi, then we obtain 

AT \I/Cor \I/Cor 

'^cord cor obs,i obs,i 

^ ~ 77 ~ N \]/raw ~ VT/raw ' V ) 

Ik ^^raw^obs,i ^obs,i 

where 'fj^^g ^ is the corrected SFR as shown in ROO. Here, we fix the data size as A'cor = A^raw = A". 
As long as Acor ~ Arawi the statistical significance is not affected by this procedure. Therefore we 
have the model value 

*i = Ci^T' , (17) 
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where 



1 rn+i 

= / ^(T)dr. (18) 



In addition, ROO note that the error bar is a Poisson fluctuation, thus the z-th error bar, crj, can be 
described as 

(„.vT/Cor \ 1/2 
. (19) 

Considering the above, we observe 

" obs,j 

k ■ 



k fvl/cor _ ^r, . 



^cor c,^''^?'" 
obs,^ -^^ '^i^obs,t 



k f\I/cpr _ vI/•^2 



E ''.2 • (20) 



• 1 

1=1 



For the numerical convenience, we fix ^'o^sj = = at r = 0, and the other 37 points are used 
for the inference (i.e., k = 37). 

As we will see later, Tgf and Tin are not determined independently. Thus, we will fix one of the 
parameters. Since Tgf is extensively investigated by Kennicutt et al. (1994), we first fix Tgf. Then, 
Tin is determined by minimizing xo^- As representative values for the gas consumption timescale, 
Tgf, we choose 0.4, 0.7 and 1.0 (tgf = 6.0, 11, and 15 Gyr, respectively)^. The best-fit Tin and xo^ for 
each Tgf are listed in Table 1. The best-fit Tin's are 1.5, 0.8, 0.6, respectively, and xo^ = 172. The 
three cases are presented in Figure 2. We see that the three parameter sets describes an almost 
identical SFH. This means that the two parameters, Tgf and Tin are strongly correlated and it is 
almost impossible to determine Tgf and Tin independently. 

We are able to test the goodness of fit between ^'obs and at the same time by evaluating 
Xo^ with respect to the x^-statistics with {k — 2) degrees of freedom. If xo^ > X^(^ ^ 2, a), then 
the hypothesis that the observed SFH is produced by the infall scenario is rejected with confidence 
level (1 — a). In fact, even if we set a = 0.01, y^{k — 2, a) = 60. Thus, we conclude that the 
data of ROO is not produced by the "classical" infall model. This clearly indicates the fact that the 
Galactic SFH is not continuous, but strongly intermittent or variable. 



^If we set Tsi < 0.4, we obtain an unreasonably large Tm(> 2) for the best fit parameter. 
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4. DISCUSSIONS 

4.1. Metallicity and G-dwarf Problem 

One of the prime motivation for infall models is that they provide a physicahy reasonable 
way of solving the G-dwarf problem (Pagel 1997 and references therein). Thus, in this section, we 
examine the chemical evolution of the Galaxy. 

Under the one-zone treatment and the instantaneous recycling approximation, the time evo- 
lution of the abundance of the heavy element, whose species is labeled by z (z = O, C, Si, Mg, Fe, 
...), is expressed as 

^ = -Xii^ + Ei + XfFit), (21) 

where Mi, Xi, Ei and X| are the total mass of heavy element i, the abundance of i (i.e., Xi = 
Mi/Mg), the total injection rate of element i from stars, and the abundance of the infalling gas, 
respectively (Tinsley 1980). Here, is expressed as 

Ei = {RXi + Yi)i;, (22) 

where R is defined in equation (2) and Yi is the mass fraction of the element i newly produced and 
ejected by stars; in other words, 

ririu 

Yi = mpi{m)(p{m) dm, (23) 

where Pi{m) is the fraction of mass converted into the element i in a star of mass m. Adopting the 
Salpeter's IMF and the stellar yield by Maeder (1992), Yq = 1.8 x 10~^, where we adopt the oxygen 
as a tracer of heavy elements (i.e., i = O; Liscnfeld &; Ferrara 1998; Hirashita 1999). Combining 
equations (1), (3), (7), (21), and (22), we obtain 

dX, ^f^_ X^:zAe-r/n._ (24) 

Since the time-evolution of /g is solved in equation (8), equation (24) can be integrated to obtain 
Xi as a function of r(= I/Tq). We hereafter choose oxygen as a tracer of the metallicity; i.e., 
i = O and assume Xl = O.OlXo,©. For a more detailed modeling, we should include the metal 
enrichment within the Galactic halo (e.g., Ikuta & Arimoto 1999). 

The result is shown in Figure 3 for the three parameter sets in Table 1, where -^o,© = 
0.013 represents the solar oxygen abundance (Whittet 1992, p. 42). The solid, dashed, and dash- 
dotted lines present the cases of (rgf, Tin) = (0.4, 1.5), (0.7, 0.8), and (1.0,0.6), respectively. We 
also present the observational data of age-metallicity relation by Rocha-Pinto et al. (2000b) (see 
their Table 3), where we assume that 0.7[Fe/H] = [0/H] = log(Xo/-^o,0) by using Edvardsson 
et al. (1993). The true parameter set which reproduce the observational data points seems to lie 
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between the solid line (rgf, Tin) = (0.4,1.5) and (rgf, Tin) = (0.7,0.8), though the discrepancy at 
the lower- metallicity side is prominent. But if the infalling gas is more enriched than have been 
assumed, this discrepancy might be resolved. We should also note the large scatter of the relation 
(Fig. 13 of Rocha-Pinto et al. 2000b) and the uncertainty in the above relation between [Fe/H] and 
[0/H]. 

The infall timescale is larger than that in Sommer-Larsen & Antonuccio-Delogu (1993), who 
gave 3.4 Gyr. The difference comes not only from the different way of modeling but also from the 
no prominent decline of the recent SFR as presented in Figure 2. 

The G-dwarf problem is also tested along with our model. The probability distribution function 
P(logXo) of the metallicity is calculated from our model as 



P(log Xo) d log Xo = * dr. 
Using equations (24) and (25), we obtain the following analytical expression for P: 



P(logXo) = (In 10) ^{t)Xo{t) 



(25) 



(26) 



where r is a function of Xq and the functional form is determined by solving equation (24). In 
comparing the distribution function with the observational data, we should take into account the 
scatter of the data. Here, we simply convolve P with a Gaussian kernel as 



/oo 
-oo 



1 



U 



27ra 



exp 



(logXo 



2a2 



dn. 



(27) 



where we adopt a = 0.1 to compare with Rocha-Pinto & Maciel (1996). We adopt these data 
because we would like to use a sample of G-dwarfs, whose lifetime is as long as the age of the 
universe. The result is shown in Figure 4. The solid, dashed, and dash-dotted lines present the 
cases of (rgf. Tin) = (0.4, 1.5), (0.7,0.8), and (1.0,0.6), respectively. The histogram shows the data 
by Rocha-Pinto &; Maciel (1996), where 0.7 [Fe/H] = [0/H] is assumed and the integrated number 
is normalized to unity. We see the same trend as Figure 3: Again we see that the true parameters 
lie between the solid line and the dashed line. As mentioned above, we should include the detailed 
enrichment process within the Galactic halo as e.g., in Ikuta & Arimoto (1999). However, since 
the aim of this paper is to examine the intermittent SFH proposed by ROO, we do not examine the 
chemical evolution further. 



4.2. Comment on the Instantaneous Recycling Treatment 

The instantaneous recycling approximation is adopted only because it provides a very con- 
venient way to obtain an analytic SFH, which is easily applied to the statistical analysis in §3.3. 
However, if we fix tgf, the instantaneous recycling treatment tends to overestimate the rate of gas 
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consumption to form stars due to the instantaneously recycled gas, which should be recycled with 
a delay in the realistic situation. Thus, we here examine the effect of the delayed recycling. We 
examine the opposite extreme case where the effect of the delayed recycling is significant: The gas 
is never recycled. In other words, we examine the case of R = 0. 

In Figure 5, we present the result for R = for the three parameters (rgf. Tin) = (0.4, 1.5), 
(0.7, 0.8), and (1.0, 0.6) with the thick lines. The thin lines represent the case where R = 0.32; 
i.e., the same as Figure 3. We adopt the same normalization for both sets of lines (equation 14). 
Comparing the thick and thin lines, we see the more decline of the SFR in the recent epoch in the 
non-recycling case than the thin lines. This is because the gas consumption timescale is shorter 
in the non-recycling case (Kennicutt et al. 1994). However, we stress that the variation of the 
observed SFH (ROO) is so large that the difference between the instantaneous and non-instantaneous 
treatments does not affect the conclusion that the SFR is violently variable. 

4.3. Two-Component Model of the Star Formation Rate 

Since the observed SFH cannot be reproduced by a simple infall model, we introduce a stochas- 
tically varying residual and modify the model as follows: 

*obs(r) = ^'infall(T)+£, (28) 

where *infaii(i) is the estimated best-fit infall model in the above discussion, and the residual, e, is 
produced by a probability distribution with zero mean and dispersion a^. We simply assume that 
is time-independent. Since e is estimated for each of the 37 bins, we can show a distribution of 
£ (Figure 6). The estimated sample dispersion = 0.22 {a = 0.47). Considering the uncertainty 
in the age estimation in ROO, the time-variation of the SFR is blurred with the uncertainty, as 
commented in ROO. Thus, the above value of a can still be underestimated. Therefore, we conclude 
that the a is at least comparable to (perhaps larger than) ^'infaii- 

If we assume that the large variation of the SFR is typical of spiral galaxies, the star formation 
activities of them should show a variety. The large value of the variance (cr) is consistent with 
previous works that suggested the variety of star formation activities of spiral galaxies (Kennicutt 
et al. 1994; Tomita ct al. 1996; Devereux & Hameed 1997). Furthermore, the kurtosis of the residual, 
K = —0.81, which means that e is distributed fiatly and is not strongly concentrated around the 
mean. Indeed, in the Figure 8 of Tomita et al. 1996, there seems to be little concentration of star 
formation activity around the mean, which is consistent with the flat distribution in Figure 6. 

4.4. Possible Mechanisms for the Variation of SFH 

Here, we consider what mechanisms are possible for the violent time-variation of SFR in the 
Galactic disk. We mention the following two mechanisms. 
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The first possibility is that the infall into the Galactic disk is a stochastic process. If an 
infalling gas is in the form of a cloud or a small-sized galaxy, such gas may induce a burst of 
star formation and increase SFR instantaneously. Indeed, infall of small galaxy seems to occur 
frequently seeing that the Sagittarius dwarf galaxy is now infalling to the Galaxy (Ibata, Gilmore, 
& Irwin 1994). The high- velocity clouds (Wakker & van Woerden 1997) may fall into the Galactic 
disk stochastically and induce stochastic bursts of star formation. 

The second possibility is that the interstellar medium is a non-linear open system. A non- 
linear open system often shows a limit-cycle behavior of physical quantities (Nicolis & Prigogine 
1977). The application of the non-linear-open-system model to the interstellar medium is described 
in Ikeuchi & Tomita (1983). According to their model, the fractional mass of the cold component 
(Xc) can oscillatory change. Recently, Kamaya &: Takeuchi (1997) suggested that the SFR varies 
oscillatory if the Schmidt law of SFR oc is assumed, where Xc is the mass ratio of the cold gas 
to the whole gas (see also Ikeuchi 1988 for a review and Hirashita & Kamaya 2000 for a recent 
development on this theme). 

5. SUMMARY 

In the context of the infall model, we re-analyzed the SFH of the Galaxy derived observationally 
by ROO. We test to examine whether variation of the star formation rate proposed by ROO is 
significant. 

We first statistically infer the timescales of the gas infall into the Galactic disk (tin) and 
that of the gas consumption to form stars (tgf)- Since each of two timescales are not determined 
independently, we first fix tgf. Then, ti^ is determined so that we minimize x^- Consequently, we 
choose three parameter sets: (igf [Gyr], ti^ [Gyr]) = (4.5, 33), (10.5, 9.0), and (15, 7.5), where we set 
the Galactic age as 15 Gyr. All of the three cases predict an almost identical SFH. The parameter 
set that seems to fit best to the age-metallicity relation and the metallicity distribution of the 
G-dwarfs is (tgf [Gyr], t-m [Gyr]) = (10.5, 9.0). The infall timescale is larger than that in Sommer- 
Larsen & Antonuccio-Delogu (1993), who gave 3.4 Gyr. The difference comes not only from the 
different way of modeling but also from the no prominent decline of the recent SFR as presented 
in Figure 2. 

Next, we test the intermittency (or violent variability) along with the smooth SFH suggested 
from the infall model. The large value of statistic supports the violent variation of the SFH. 
Then, we interpret the observed SFH with the two components; underlying smooth component 
described by the infall model, and violently variable component. We find that the variation of the 
latter component is comparable to the former. 

As a test of the models, we also examine the age-metallicity relation and the metallicity dis- 
tribution of the Galactic stars. Consequently, we observe that the degeneracy of the three param- 
eters are resolved and find that the true parameter set seems to lie between (tgf [Gyr], tin [Gyr]) = 
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(10.5, 9.0) and (4.5,33). 

Finally, two physical mechanisms for the variable SFH are suggested. One is the stochastic 
infall of the clouds or small galaxies, and the other is the non-linear oscillation of the SFR due to 
the limit-cycle behavior of the fractional mass of the cold-gas component. 
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A. SMOOTHING WITH MOVING MEDIAN 

In this appendix, we explain the details of the EDA smoothing based on the moving median. 
The procedure is as follows: 

1. We take the moving median of three sequential values of the original data. For the boundary 
data, we use the median of the following: 

xi = median(x2 — 2(a;2 — xa), xi, X2) , (-^1) 
— median(xn + 2(x^_i — Xn—2) 1 Xn—i-, Xn 

) , (A2) 

where Xi is the smoothed value. 

2. We perform the splitting of the plateau of the median-smoothed data sequence, i.e. when we 
find the same value Xi and Xj+i, we substitute 

Xi = medisin{xi-i + 2{xi-i- Xi-2),Xi,Xi-i), (A3) 
Xi+i = median(a;i+2 - 2(a;i+3 - a;i+2),Xi+i,Xi+2). (A4) 

into Xi and Xi+i, respectively. 

3. We, then, take the mean of Xi-i and (we denote them as Xj). Finally, we derive the 
mean of Xi and Xj . 

This method is often used in the econometric and biometric researches, and provides successful 
results. 
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Table 1: Examined Parameters and xo^- 





tsf (Gyr) 


'in 


tin (Gyr) 


Xo 


0.4 


6.0 


1.5 


23 


172 


0.7 


11 


0.8 


12 


172 


1.0 


15 


0.6 


9 


172 



"The gas consumption timescale {tst) normalized by Tq = 15 Gyr. 
''The gas infall timescale (tin) normalized by Tq = 15 Gyr. 
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Fig. 1. — The smoothed star formation history of the Galaxy. In the upper panel, we show 
the smoothing result by moving average. The dot-dot-dot-dashed line is the smoothed SFR with 
smoothing kernel width of 1.2 Gyr, strong dashed line is the smoothed SFR with 2.0 Gyr kernel, 
and dot-dashed line is the smoothed SFR with 2.8 Gyr kernel. In the lower panel, the dashed line 
represents the moving average of 2.0 Gyr kernel, and the solid line is the smoothed SFH by the 
moving median method of the exploratory data analysis. We observe that the 1.2-Gyr kernel is not 
sufficient to smooth the varying SFH. 



-17- 



2.5 




0.0 I . . . . I . . . . I . . . . I 

5 10 15 

To - t [Gyr] 



Fig. 2. — Star formation histories for the parameter sets in Table 1. The dotted hne is the 
observed star formation history in ROO. The sohd, dashed, and dash-dotted hnes present the cases 
of (rgf, Tin) = (0.4, 1.5), (0.7, 0.8), and (1.0, 0.6), respectively, though it is difficult to distinguish 
the three lines. 



0.5 




Fig. 3. — Time evolution of oxygen abundance Xq for the three sets of the parameters in Table 
1. The solid, dashed, and dash-dotted lines present the cases of (rgf, rin) = (0.4, 1.5), (0.7, 0.8), 
and (1.0, 0.6), respectively. Contrary to Fig. 2, we can easily distinguish the three lines. The true 
parameter set seems to lie between (rgf. Tin) = (0.4,1.5) (the solid line) and (rgf. Tin) = (0.7,0.8) 
(the dashed line), though uncertainty and discrepancy exist. 
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Fig. 4. — Distribution of the G-dwarf metallicity. for the three sets of the parameters 
in Table 1. The soHd, dashed, and dash-dotted hnes present the cases of (rgf, Tin) = 
(0.4, 1.5), (0.7, 0.8), and (1.0,0.6), respectively. The histogram shows the data by Rocha-Pinto 
& Maciel (1996). Again the true parameter set seems to lie between (rgf, Tin) = (0.4, 1.5) (the solid 
line) and (rgf. Tin) = (0.7,0.8) (the dashed line), though uncertainty and discrepancy exist. 




Fig. 5. — Star Formation History in the non-recycling model {R = 0) for the three parameters {thick 
lines) as an extreme case, where the timcscalc of the recycling is very large. The star formation 
history in the instantaneous recycling model is also shown {thin lines). The solid, dashed, and 
dash-dotted lines present the cases of (rgf. Tin) = (0.4, 1.5), (0.7,0.8), and (1.0,0.6) , respectively. 
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Fig. 6. — Distribution of the residual component of the star formation rate, e. 



